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Gravitational waves from coalescing compact binaries are one of the most promising sources for 
detectors such as LIGO, Virgo and GEO600. If the components of the binary posess significant 
angular momentum (spin), as is likely to be the case if one component is a black hole, spin-induced 
precession of a binary's orbital plane causes modulation of the gravitational-wave amplitude and 
phase. If the templates used in a matched-filter search do not accurately model these effects then 
the sensitivity, and hence the detection rate, will be reduced. We investigate the ability of several 
search pipelines to detect gravitational waves from compact binaries with spin. We use the post- 
Newtonian approximation to model the inspiral phase of the signal and construct two new template 
banks using the phenomenological waveforms of Buonanno, Chen and Vallisneri [A. Buonanno, 
Y. Chen and M. Vallisneri, Phys. Rev. D67, 104025 (2003)]. We compare the performance of these 
template banks to that of banks constructed using the stationary phase approximation to the non- 
spinning post-Newtonian inspiral waveform currently used by LIGO and Virgo in the search for 
compact binary coalescence. We find that, at the same false alarm rate, a search pipeline using 
phenomenological templates is no more effective than a pipeline which uses non-spinning templates. 
We recommend the continued use of the non-spinning stationary phase template bank until the false 
alarm rate associated with templates which include spin effects can be substantially reduced. 

PACS numbers: 04.30.Db,04.25.Nx,04.80.Nn,95.55.Ym 



I. INTRODUCTION 



In 2005-2007, the Laser Interferometer Gravitational- 
Wave Observatory (LIGO) recorded two years of data 
at design sensitivity [l| and the LIGO, Virgo 01 and 
GEO600 0] detectors now form a world-wide network of 
broad-band gravitational-wave observatories. The LIGO 
and Virgo detectors are scheduled to resume operations 
in summer 2009 with a factor of ~ 2-3 sensitivity increase 
over previous observations. The gravitational waves 
emitted during the inspiral and merger of binaries con- 
taining neutron stars (NS) and/or black holes (BH) are 
a primary target of this network. Binary neutron stars 
(BNS) can be observed up to 35 Mpc (70 Mpc) in the 
Initial (Enhanced) LIGO detectors and up to 450 Mpc in 
the Advanced LIGO detectors, which will begin observa- 
tions in 2015 Binary black holes (BBH) with lOA/© 
components should be visible at 170 Mpc (350 Mpc) in 
the Initial (Enhanced) LIGO detectors, increasing to 2 
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Gpc in Advanced LIGO [4|. Population synthesis cal- 
culations constrained by radio observations of BNS sys- 
tems containing pulsars predict BNS detection rates be- 
tween 10~^-1 yr~^ for Enhanced LIGO and 0.4-400 yr~^ 
for Advanced LIGO, with the most likely values being 
0.1 yr~^ and 40 yr^^, respectively [1, H, B]. Much less 
is known about the detection rates of BBH and NS-BH 
coalescences, although it is plausible that Enhanced (Ad- 
vanced) LIGO rates could be as high as 1 (300) yr""'^ for 
NS-BH binaries and 20 (4000) yr"! for BBH [4|,0,[3- 

The sensitivities listed in the preceding paragraph are 
optimal: they assume accurate knowledge of the signal 
waveform in order to construct matched filters which can 
extract gravitational-wave signals buried in the noisy de- 
tector data |T3| ■ The gravitational waveform from the 
inspiral of two compact objects has been calculated us- 
ing the post-Newtonian (PN) approximation, which uses 
the characteristic velocity of the binary (f/c ) as an ex- 
pansion parameter 0, M, [H, [H, [H, M, [13, 111, E, M, 

EH, [13, [11, [13, [111 • Ongoing comparisons of PN wave- 
forms with numerical simulations of binary black holes 
have thus far confirmed the accuracy of the PN solution 
in the late stages of inspiral [H, [13, [111, although optimal 
searches for sources with total mass > 30 Mq in the first- 
ij^imeiEalij on detectors require waveforms that also model 
the merger and ringdown (29l . [30j . If the components 



of the binary have negligible instrinsic angular momen- 
tum (spin) then it is straightforward to construct a bank 
of matched filters, parameterized by the two component 
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masses of the binary, and use these filters to search for 
signals [Sll, [H, [H, [SJ] . However, if the components of a 
binary are spinning, then these spins can couple with the 
orbital angular momentum of the binary and with each 
other to cause amplitude and phase modulation of the 
gravitational waveform [35l | . Attempting to detect gravi- 
tational waves from spinning binaries with non-spinning 
templates will result in a sub-optimal search and a corre- 
sponding reduction in the detection rate [H, [s^l ■ Since 
it is possible that a large fraction of astrophysical black 
holes have considerable spin (ssl . [s^ , it is important to 
consider the effect of spin in searches for gravitational 
waves from BBH and NS-BH coalescences. However, op- 
timal filters for spinning binaries are characterized by a 
much larger number of parameters than the ones for non- 
spinning binaries, complicating placement of filters in the 
bank and considerably increasing the computational cost 
of searches. 

To mitigate the computational problem without com- 
promising the sensitivity of the search, a phenomenologi- 
cal family of templates was proposed by Buonanno, Chen 
and Vallisncri (we refer to these templates as BCV- 
spin). Filters constructed from these templates are de- 
scribed by only four parameters and have good overlaps 
with the full PN waveforms Moreover, constructing 
a bank of filters using BCV-spin waveforms is straight- 
forward, if cumbersome [4l|, [43| . The first searches for bi- 
nary black hole signals in LIGO data used non-spinning 
templates [l^E^], however BCV-spin templates were re- 
cently used to search for BBH and NS-BH signals with 
spin in data from the third LIGO science run [i^]. The 
sensitivity of the search described in [i^l was not as good 
as the main results of Ref. [i^ might suggest. This 
was primarily due to the the response of the BCV-spin 
template to the non-gaussian noise transients present in 
real gravitational-wave detector data and the increase in 
the number of degrees of freedom associated with the 
detection statistic (due to the larger search parameter 
space) [42| . This was already anticipated in [40| ; here we 
provide a quantitative analysis. 

In this paper, we present an improvement to the search 
pipeline described in [4^ . by constructing banks that 
are much better suited to the BCV-spin template family. 
We compare the sensitivity of this search to the search 
for gravitational waves from compact binary coalescence 
with non-spinning templates in LIGO data [1^, [s^l • Our 
main conclusion is that while the BCV-spin templates 
have rather good overlaps with the target waveforms, the 
current search pipeline needs further improvements be- 
fore any gains from these increased overlaps can be re- 
alized. The false alarm rate of BCV-spin filters in real 
detector data is larger than that of a non-spinning search. 
This makes a search using BCV-spin templates less sen- 
sitive than a non-spinning search when looking for bina- 
ries with spin, since one has to use a higher detection 
threshold to obtain the same false alarm rate. The re- 
sults of this work were used to guide the decision not to 
implement the BCV-spin search on data from the fifth 



LIGO science run and instead to use non-spinning fil- 
ters to search for binaries with spin (45j . The motiva- 
tion for this decision was summarized in an appendix of 
Ref. [iH and this paper can be seen as a companion to 
that work. Here we present a detailed account of how the 
BCV-spin banks were constructed, and how the compar- 
isons between the BCV-spin and non-spinning searches 
were performed. 

This paper is organized as follows. In Sec. |TT] we 
give a description of our target signals, which are post- 
Newtonian waveform models for signals from spinning 
black hole binaries, followed by a summary, in Sec. Illli 
of the phenomenological BCV-spin templates of Ref. [4C| . 
In Sec. IIVI we review the construction of the template 
bank used in Ref. [4^ and present two new methods 
to construct BCV-spin template banks, relaxing the as- 
sumptions used in Ref. [i^ : a "square-hexagonal" place- 
ment which generalizes the hexagonal placement devel- 
oped in (46j to a higher-dimensional template manifold, 
and a stochastic placement technique proposed in [i^] . In 
Sec. |V] we compute matches of these banks with target 
waveforms and compare the results with those obtained 
from two-dimensional template banks based on the sta- 
tionary phase approximation (SPA) [H, [4^ and hexag- 
onal tem plat e placement used in LIGO's non-spinning 
searches lH Sj, El, [Hill IH, HI- We compare the 
detection efficiency of the spinning banks with that of 
SPA banks in Sec. IVII followed by concluding remarks in 
Sec. I VIII Throughout this paper, we set G ~ c = 1 unless 
otherwise stated. 



II. POST-NEWTONIAN WAVEFORMS FROM 
SPINNING BINARIES 

Depending on their birth spins, BHs in binaries could 
accumulate significant spin through accretion [H, [s^ . 
There is much uncertainty concerning the equation of 
state of a neutron star, but most models place an upper 
limit J/Ml< 0.7 on the spin, above which the star would 
break up |55| . There is also an upper limit for the spin 
of a black hole due to torque caused by radiation from 
the accretion disk getting swallowed by the BH, leading 
to an expected bound of J/AP < 0.998 Most of 

the modeling of spin evolution in compact binaries has 
been confined to NS-BH systems, in which case the spin 
tilt with respect to the orbital angular momentum can be 
considerable [H,!!^; this may also be the case for BBHs. 

PN theory has achieved great success in modeling 
the adiabatic, quasi-circular phase of inspiral, during 
which the fractional change in the orbital frequency 
over each orbital period will be negligible (see, e.g., 
Ref. [13 for a review). The orbital phasing has been 
calculated 
tion) 

while the gravitational-wave amplitude for non-spinning 
binaries has been calculated to order (v/c)^ [Tsl. [l9l. [20| . 
The effect of spin on the gravitational-wave phasing is 
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known to order (■y/c]f, [2TI, [H, [H, [ll] and to order [v/c)^ 
for the amplitude [2^. However, since the matched fil- 
ter is most sensitive to the phase evolution of the bi- 
nary, template waveforms amplitudes are typically com- 
puted only at leading order in amplitude (the restricted 
waveform). Spin-orbit interaction enters the phasing at 
1.5PN and 2PN order and spin-spin interaction at 2PN 
order. Spin effects influence the evolution of the orbital 
frequency as a function of time. Including these effects, 
the adiabatic evolution of the orbital frequency oj{t) is 
given by^ 
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where Ljv is a unit vector in the direction of orbital angu- 
lar momentum (and hence the unit normal to the orbital 
plane of the binary), Si^2 are the spins, xi, 2 = |Si,2|/?7ii,2 
with 7711^2 the component masses, 81^2 = Si_2/|Si^2|, 
M = 777,1 + m-2 is the total mass, and 77 = 17111112/ M'^ 
the symmetric mass ratio. 7 = 0.577 ... is the Euler- 
Mascheroni constant. The spins and the direction of the 



^ At the time this work was started the spin-orbit term at 
2.5PN I23I I24II was not yet known and so is not included here. 



orbital angular momentum evolve according to [22l . |3E 
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where S = Si + S2. The dynamics of the binary is 
governed by the nonlinear, coupled differential equations 
([T])-(|31)- It will not be possible to solve these exactly, but 
they can easily be treated numerically. 

By numerically evolving uj{t) one can obtain the orbital 
phase. 



idt. 



(5) 



which can be substituted into the usual expressions for 
the restricted FN waveform polarizations [6J|. In the 
case of spinning binaries, we need to take into account 
the time-dependence of the amplitudes through the in- 
clination of the orbit with respect to the observer. The 
plus and cross polarizations of the gravitational wave are 
given by 

h+{t) = -[l + (L-n)2j cos(2$(i)), 

hy,{t) = -2(L-n) sin(2$(i)), (6) 

where the unit vector n points from source to detector. 
The detector strain is 



h{t) ^ F+{t)h+{t) + F^{t)h^{t), 



(7) 



where F+^^{t) = F+^x (a, 'l'(i)) are detector antenna 
factors which depend on the the right ascension and dec- 
lination of the source and a time-dependent polarization 
angle *(t) (see, e.g., Ref. [ssj). 

As suggested by Eq. (|4]), the direction of orbital an- 
gular momentum and hence the plane of the inspiral will 
undergo precession, the effect being more pronounced for 
asymmetric systems. It will also be more prominent if the 
spins are large, and if they are significantly misaligned 
with the orbital angular momentum. The time evolution 
of spins and angular momentum will affect the phasing of 
the waveform through Eqns. (H])-®, and the precession 
of the orbital plane will modulate the amplitudes of the 
wave polarizations in Eq. ([5]). The waveforms given by 
Eqns. ([T])-([71) will be the "target-signal waveforms" for 
testing our template banks. 
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III. A DETECTION TEMPLATE FAMILY FOR 
SPINNING BLACK HOLE BINARIES 

The frequency-domain phenomenological detection 
template family proposed in Ref. [i^l is designed to cap- 
ture spin-induced amplitude and frequency modulation 
in an approximate way. Specifically, for gravitational- 
wave frequencies / > 0, the BCV-spin template is 



e(/cut-/) 



3 

E 



Here to is the time of arrival, Q(x) is the usual Hcaviside 
step function and /cut is an upper cut-off frequency be- 
yond which the waveform is unlikely to be close to a true 
signal (due to breakdown of the adiabatic approximation 
to the inspiral regime). The detection statistic will be 
maximized analytically over the parameters a'l, . . . , in 
the linear combination ([8|), as well as over to', these pa- 
rameters are referred to as extrinsic parameters because 
they do not need to be explicitly searched over. 

The waveforms hj(f), j = 1, . . . , 3 are basis templates, 
which take the form 



i*NM(/) 



(9) 
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and P captures the effect of spin-induced amplitude mod- 
ulation. The (non- modulated) phase <&nm(/) takes the 
form^ 



*NM(/) = .r'/'(^0 + ^'3/). 



(11) 



It will not be possible to analytically maximize the de- 
tection statistic over the parameters ipo, ips, and /?, and 
these must be explicitly searched over using a bank of 
templates; they are referred to as intrinsic parameters. 

It will often be useful to approximately identify the 
intrinsic parameters with the physical masses and spins 
of a compact binary. By rela ting tpo and ip3 to the OPN 
and 1.5PN phase coefficients [63], one has the correspon- 
dences 

^0 ^ ^(7rM)-5/3, ^3 ^ 



(12) 



Similarly, the parameter P can be related to the rate of 
precession by [35[ 
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^ What is called f/'s here was denoted ipi in [4^ 



We stress that these mappings are only approximate, and 
for a given physical signal, the detection template that 
matches best may correspond to values of {ipo, ips, /3) that 
differ significantly from the ones suggested by the iden- 
tifications above. 

The identifications allow us to make a choice for 
/cut- In the limit where one component mass goes to 
zero while total mass M remains fixed, and assuming 
zero spins, the frequency of last stable orbit (LSO) of 
a test mass in the Schwarzschild spacetime is given by 
/lso(M) = (63/27rM)-i. For simplicity we set /cut = 
fhsoiM), where M = —ip^/^lGTr^ipo) is computed from 
the correspondence (fT2|) . 

Next one constructs an orthonormal basis from the ba- 
sis templates ([9]) with respect to the usual inner product 
for waveforms a, b on the template manifold given by 



(a, b) = m 



/s 



a{f)b*{f) 

Snif) 



df, 



(14) 



where tilde denotes a quantity computed directly in the 
frequency domain (as in the case of the BCV-spin tem- 
plates) or the Fourier transform of a time-domain quan- 
tity (such as the waveforms h{t) given in Eq. ([7])). Sn{f) 
is the one-sided power spectral density (PSD) of the de- 
tector data, and fs is some lower cut-off frequency asso- 
ciated with the detector; in the case of initial LIGO one 
sets fs = 40 Hz. The orthonormalization of the basis tem- 
plates can be effected using the Gram-Schmidt procedure 
as in [l^j. In addition one demands that the templates 
themselves are normalized (denoted by h): {h,h) ~ 1. 
This leads to the requirement 
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(15) 



where the a^, j = 1, . . . , 6, are the coefficients of h when 
expressed into the orthonormal basis of templates result- 
ing from the Gram-Schmidt procedure. 

Finally, the signal-to- noise ratio (SNR), which is used 
as the BCV-spin detection statistic, is given by 



ll/2 



max(s, h[to, ctj]) 



(16) 



where s represents the detector data stream, and the 
maximization over the aj is subject to the constraint 



IV. TEMPLATE BANKS FOR SPINNING 
BINARIES 



The template waveforms h may not exactly model 
gravitational- wave signals s. The loss in SNR due to 
differences between the template and signal waveforms is 
quantified by the fitting factor J- [ssj . If s is a signal 
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waveform and h a template waveform, then 

= max(s, /i), (17) 

h 

where hat denotes normaUzatfon: (s,s) = {h,h) — 1. 
1 — is the fractional loss in SNR resulting from the 
use of sub-optimal template waveforms rather than the 
true signal waveforms. Since we do not a priori know the 
intrinsic parameters of any gravitational-wave signals we 
may detect, we decide on a target signal space and con- 
struct a discrete bank of templates to search for signals 
in this space. If hb is a normalized template waveform in 
the discrete bank and /i is a normalized waveform from 
the space used to construct the bank, then the minimum 
match M of the bank is defined to be [3l| 

M = min I max {h, hb) j (18) 

h \/ib£bank / 

A typical choice for the minimum match in gravitational 
wave searches is = 0.97. When measuring the per- 
formance of a template bank, we are interested in the 
effective fitting factor f given by [1^ 

T = _ max (s./if,) (19) 

/if,£bank 

If the signal waveforms are identical to those used to 
construct the bank, then the effective fitting factor will 
be bounded below by the minimum match. In practice, 
the true gravitational-wave signals will differ from the 
templates used to construct the bank, so the effective 
fitting factor may be smaller than the minimum match. 
The larger the effective fitting factor, the better the bank 
is at capturing the target signals. 

If the parameters A between two (normalized) tem- 
plates differ by a small amount AA, the loss in SNR can 
be related to a distance defined by a metric gij given by 

(J/i(A),/i(A-F AA)^ ~ l-5yAA*AA^ (20) 

where 

_\l dh_ dh_\ 

The standard method of constructing a template bank 
then consists of computing this metric in the intrinsic 
parameter space of a waveform family and using it to 
place templates such that the distance between any tem- 
plate waveform and the nearest template in the bank is 
greater than the desired minimum match M . In searches 
for non-spinning binaries, the intrinsic parameters of the 
templates are the just component masses (mi, 7712) of 
the binary. In practice, we re-parameterize the wave- 
forms using the chirp times (ro(mi, 7712), r3(mi, TO2)) 
With respect to these variables the metric is almost Eu- 
clidean, and so template placement using the metric 



gij becomes a straightforward two-dimensional hexago- 
nal packing problem [ssj . 

As described in Sec. [H a search for gravitational waves 
using the BCV-spin templates has been performed in S3 
LIGO data The metric used in that analysis was 

computed using the "strong modulation approximation" 
where one assumes that the binary precesses many times 
while emitting in the most sensitive part of the detector's 
band. This allows one to treat the basis templates of 
Eq. ([9]) as orthonormal, simplifying the calculation of the 
metric. However, the resulting template banks were only 
appropriate for fairly low-mass, asymmetric systems. We 
now present an improved algorithm for constructing a 
metric in which the assumptions of [4^ are dropped. In 
our case, the parameters of the waveform are 

A= (to,ai,...,a6,V'o,V'3,/3)- (22) 

The detection statistic can be maximized over the ex- 
trinsic parameters tg and ai, . . . , ag, which, as shown in 
[67t . leads to a projected metric gfj°^ which only mea- 
sures distances in the (V'Oj V'S: Z^) directions. However, 
the components of 5^^^°^ will still depend on the aj . This 
residual dependence on extrinsic parameters can be re- 
moved as follows: 

1. Introduce some fiducial distance Asq; 

2. Specify a large number of unit vectors (in the co- 
ordinate sense) n in (t/jq, ?/'3, /3) space; 

3. For each n, numerically maximize the metric length 
Asf, computed from gfj°\ over values of the aj 
consistent with the constraint psp : i.e., 

A4= max gf'°''{a^,P)h'h^; (23) 

4. Rescale each vector h by defining a new vector u = 
(Aso/Asft)n; 

5. Fit an ellipsoid in parameter space to the vectors 
u; 

6. Define an "effective" metric g^^ by requiring that 
any point on the ellipsoid is at effective metric dis- 
tance Aso from the template we started with. 

Note that this construction is independent of the fidu- 
cial length scale Aso- In what follows, is the metric 
we will use to satisfy the criterion (|18p through the rela- 
tionship (PO]) . An property of g^j^ is that it is essentially 
independent of ipQ and "03 and only has a weak depen- 
dence on (3. 

It is important to note that given a short straight line 
segment in coordinate space with coordinate length AA, 
by construction g^^ associates almost the largest possi- 
ble metric length to it consistent with the family of met- 
rics gfj°'' (ctj , /?) parametrized by the aj . When generat- 
ing template banks, in practice one specifies a minimum 
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match which will then be used together with the metric 
to determine the spacing of templates. Since g'^^ is too 
conservative in assigning lengths, neighboring templates 
will tend to have a larger match than needed, and the 
true minimum match defined by (|18p will always be sig- 
nificantly larger than what was originally intended. As 
we shall see below, setting an a priori value of = 0.8 
will be more than enough for a bank to obtain high over- 
laps (> 0.9) with target waveforms. 

We would like to capture signals from binaries whose 
component masses are in the interval [1, 35]Mq, with to- 
tal masses M < 35 Mq. We do not need to worry about 
capturing BNS signals, since spin does not have a signif- 
icant effect on waveforms from those sources. However, 
our template bank should have good overlap with NS- 
BH and BBH signals. Taking neutron star masses to lie 
between 1 Mq and 3 il/0 and black hole masses to be 
larger than 3Af0, we impose M > AMq. To capture 
these signals, we want an appropriately chosen bound- 
ing box in [tpQ,^^, (3) within which to place templates. 
Such a box can be specified using the correspondences 
P^ - (fT^ . The suggested intervals for {tpo^i's) are then 
roughly 

^/-o G [8 X 10^5 X 10^]Hz^/^ 

lAa e h3 X 10^ 10] Hz2/^ (24) 

where the upper bound for tp^ has been chosen gener- 
ously. As to /3, the correspondence suggests that 
P < 150 Hz^/^ should suffice, but to have good matches 
with a variety of physical signals, here too it turned out 
to be better to have a larger upper bound: 

/3 e [1,4 X 102]Hz^/^ (25) 

We now present two methods for constructing template 
banks for BCV-spin templates which cover this space. 

A. Square-hexagonal template bank 

The metric g°j^ depends only on /3, so it is natural 
to first define layers of constant /3, with a spacing de- 
termined by the minimum match. Within each of the 
two-dimensional layers one can then lay out templates in 
a hexagonal pattern, which is the optimal placement in 
two dimensions. We will refer to this kind of placement as 
square-hexagonal. The construction of this bank is analo- 
gous to that described in Ref. [3 which was used to con- 
struct template banks for search for binary black holes 
in data from the third and fourth LIGO science runs [13] 
using non-spinning phenomenological templates [68| . For 
the BCV-spin templates, we have a 3-metric, which in 
each (3 layer is diagonalized by going to a new set of co- 
ordinates (i/'g,-03,/?'), where = /3. After that a hexag- 
onal placement in (tp^jip'^) can be performed as in [46| . 
As explained above, the metric is overly conservative in 
specifying distances between templates, and setting an a 
priori minimum match of A4 = 0.8 will suffice to obtain 
high matches with target waveforms. 



B. Stochastic template bank 

We now consider a different bank placement for BCV- 
spin, which we hope will reduce the overcoverage of the 
parameter space that is unavoidable with the square- 
hexagonal placement method defined above. This will 
lead to a smaller number of templates but will yield the 
same or better matches with target waveforms, and sim- 
ilar efficiencies. This template bank is created by the 
placement of a large number of randomly distributed 
templates, followed by a "pruning" stage in which un- 
necessary t emp lates are discarded. This method is de- 
scribed in [47j and summarized below. Other, similar 
methods for creating stochastic template banks were pro- 
posed in [i^ and [ZO]. 

The stochastic placement algorithm we wish to use for 
BCV-spin is very simple. We begin by generating a very 
large number of points in the parameter space, far more 
than would be needed to fill the space. We then itera- 
tively cycle through these points, retaining a point only 
if it is not closer than some predefined metric distance 
A to the points retained in previous iterations. The re- 
maining points form our stochastically generated bank. 
Tests have shown [13] that one should begin with at least 
N^-^ points, where N would be the number of templates 
remaining after filtering, to have a good coverage of the 
parameter space after "pruning" . 

In testing this algorithm against lattice placement al- 
gorithms it was found [13] that in a 2-dimensional Carte- 
sian space the stochastic algorithm produced a template 
bank with 1.5 times the number of templates that a 
square lattice algorithm would have generated. How- 
ever, in the case of a 2-dimensional non-spinning (non- 
Cartesian) SPA bank (as described above) the stochas- 
tic algorithm was found to place ~ 10% less templates 
than the square lattice algorithm and only ~ 25% more 
templates than the hexagonal lattice placement, while 
achieving a similar degree of coverage. We emphasize 
here that this stochastic placement algorithm would be 
of most use in parameter spaces with more than 2 di- 
mensions, where lattice placement becomes significantly 
sub-obtimal. 

For the specific case of BCV-spin the templates are 
sprinkled randomly over a rectangular box in (t/jq, "031 /5) 
space using the same bounding box as in the previous 
subsection. An estimate for the number of templates that 
will be needed is provided by the invariant volume of the 
box, divided by the volume taken up by an individual 
template: 

Once again it will suffice to set an a priori minimum 
match = 0.8 (i.e., setting the A defined above to 0.2). 
Given the box in parameter space specified by (|24l25p . 
the number of sprinkled templates should then be about 
500,000. When using a larger number of initial templates. 



7 



we find that the final number of templates after pruning 
does not change significantly. With the Initial LIGO de- 
sign PSD, the number of templates for stochastic BCV- 
spin banks with M = 0.8 is about 8000; SPA banks 
with M = 0.95 have ~ 12,000 templates, and for BCV- 
spin with square-hexagonal placement and jV( = 0.8 more 
than 16,000 templates are obtained (see Table H]). 



V. COMPARISON OF BCV-SPIN BANKS WITH 
SPINNING PN SIGNALS 

We now study the performance of our banks against 
the target waveforms of Sec. |lll In particular, for a va- 
riety of target waveforms s corresponding to different 
masses and initial spins, we compute the effective fitting 
factor T of the bank for the target waveforms, as given 
by Eq. (HH). 

Fig.[l]compares the effective fitting factor of templates 
in square-hexagonal and stochastic BCV-spin banks with 
those of a non-spinning SPA bank. There is no discern- 
able difference between square-hexagonal and stochastic 
placements, but both differ significantly from the non- 
spinning SPA bank. As one would expect, the difference 
is largest for binaries with large mass ratios, although 
there is improvement also for a variety of other target 
waveform masses. Depending on masses and spins, for 
the same target waveform s, the difference in can be 
more than 25%. The medians and means for the ef- 
fective fitting factors are summarized in Table IH We 
find that BCV-spin with stochastic bank placement has 
marginally better effective fitting factors than the square- 
hexagonal bank, and both BCV-spin banks have notica- 
bly higher effective fitting factors than the non-spinning 
bank. Given the small difference between the stochas- 
tic and square-hexagonal BCV-spin banks we will subse- 
quently only consider differences between the stochastic 
BCV-spin bank and the non-spinning SPA bank. 



VI. SEARCH PERFORMANCE OF BCV-SPIN 
TEMPLATE BANKS 

The effective fitting factor of a target waveform over 
a template bank as defined in Eq. indicates how 
similar the templates are to physical signals, but when 
searching for gravitational-wave signals in real detector 
data, other factors also come into play. The effective 
fitting factor of a template bank gives us a measure of 
how the signal-to-noise ratio is reduced by not filtering 
with the true signal waveform, but to detect a signal we 
must be able to distingusih it from background noise in 
the detector. To determine the overall performance of a 
template bank, we have to consider both the effective fit- 
ting factor and the false alarm rate of the bank, i.e., the 
response of the filters to noise (both Gaussian and tran- 
sient) in the detector. Once we establish the false alarm 
rate of a search, we measure the performance of a bank by 



BCVSpin Square-Hexagonal 




mi (Mo) 



Figure 1: Effective fitting factors of 1124 target waveforms 
with templates in a square-hexagonal BCV-spin bank (top), 
a stochastic BCV-spin bank (middle), and an SPA bank us- 
ing the Initial LIGO design PSD. The black crosses indicate 
the component masses for the target waveforms; spins have 
random orientations, and 0.7 < |Si,2/nii 2I < 1- The color 
coding gives the effective fitting factor of a target waveform 
over the template bank. There is no discernable difference be- 
tween the performances of the square-hexagonal and stochas- 
tic BCV-spin banks, but both do notably better than SPA. 
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Table I: Number of templates for the SPA, square-hexagonal, and BCV-spin banks with Initial LIGO design PSD, and median 
and mean effective fitting factors T of the banks with target waveforms as in Fig. [T] 



Template 


Bank-placement 


SNR threshold 


Minimum match M 


Number of templates 


median 


rncan 


SPA 


SPA 


5.5 


0.95 


11832 


0.89 


0.86 


BCV-spin 


Square-hexagonal 


8 


0.8 


16431 


0.96 


0.92 


BCV-spin 


Stochastic 


8 


0.8 


7913 


0.96 


0.93 



its efficiency^ i.e., the bank's ability to find simulated tar- 
get waveforms injected in the noise at a given false alarm 
rate. We will establish the false alarm rates and efRcien- 
cies of BCV-spin and SPA banks by means of the data- 
analysis pipeline used in searches by the LIGO Scien- 
tific Collaboration (LSC) for inspiral signals [13, 0, 45 1, 
which is available in the LSC Algorithm Library 71 1. 
More details on this pipeline can be found in Ref. [72t. 



A. False alarm rates 

The BCV-spin detection statistic (|16p involves maxi- 
mization over six parameters (ip, and the with the 
constraint '^jCtj = 1), to be compared with only two 
for SPA. It should also be noted that the BCV-spin de- 
tection template family consists of waveforms that are 
only approximate. As we shall see below, the larger 
number of degrees of freedom will make the BCV-spin 
banks more prone to detecting instrumental noise tran- 
sients with high SNR. Both for SPA and BCV-spin one 
needs to set an SNR threshold below which no candidate 
events are accepted, and the higher false alarm rate with 
BCV-spin will necessitate setting a higher threshold. 

The pipeline used to search for gravitational-wave sig- 
nals in the LIGO detectors demands that candidate 
events be coincident in two or more detectors [t^- If 
the noise sources in our detectors are uncorrelated (as 
in the case of the two 4 km LIGO detectors), we can 
measure the false alarm (or background) rate of this 
pipeline by time-shifting the detector data by more than 
the gravitational-wave travel time between the detectors 
(11 ms) and looking for coincident triggers; such triggers 
will be due to accidental coincidence of noise alone. We 
can repeat this with time steps of, say, 10 s, and count the 
number of coincident triggers in each of the time-shifts 
to obtain a good estimate of the false alarm rate. 

Before triggers are compared between detectors, they 
are clustered together, keeping only the trigger with the 
loudest SNR within a certain time window (in our case 
4000 ms). Next, various methods can be used for declar- 
ing two clustered triggers to be coincident across detec- 
tors. Usually one demands not only coincidence in time, 
but also that the parameters of the template that gave 
the loudest SNR be similar in the different detectors. 
The simplest way of implementing this is the so-called 
box-coincidence method, whereby two triggers are con- 
sidered coincident if they occured within a certain time 
from each other (say, 100 ms), and the associated tem- 
plates have parameters that differ only within certain 



tolerances [T^l- In the case of BCV-spin, these were cho- 
sen as A^o = 40000 Hz^/^ and = 600Hz^/^ with 
no restrictions on differences in /3 [42| . 

More recently, a more sophisticated technique was de- 
veloped which has the potential to dramatically reduce 
the false alarm rate [t^ . In this method, the covariances 
between the signal parameters are used to define an er- 
ror ellipsoid in parameter space around the triggers, and 
triggers in different detectors are considered coincident 
if their associated ellipsoids overlap. In the case of SPA 
banks, the size of the ellipsoids will depend strongly on 
the region of parameter space the triggers occur in. Gen- 
erally they will be smaller for triggers associated with 
smaller masses, as waveforms will then spend more time 
in the detector band and errors will be smaller. This 
leads to a dramatic reduction in the number of spuri- 
ous coincident triggers. By contrast, the box-coincidence 
method described above uses the same parameter win- 
dows anywhere in parameter space. 

The ellipsoid coincidence method has been successfully 
implemented for SPA banks. The technique is well-suited 
for banks where the templates are simplified versions of 
target waveforms, so that one can assume template wave- 
forms to be reasonably close to physical signals. It would 
be possible in principle to implement such a method also 
for the (phcnomcnological) BCV-spin banks. However, 
in this case the metric gf^^ is basically independent of ■00 
and -03, the parameters that are most closely related to 
the masses. Hence, for a given value of /?, the associated 
ellipsoids would not differ in size across ('i/'Oi'/'a) space, 
and only their orientations would differ with p. This 
way, no great improvements can be expected in terms of 
reducing the false alarm rate. 

Tabic, [ni shows the average number and variance of 
coincident triggers between the 4 km LIGO Hanford and 
Livingston detectors for time-shifts within ~ 9 days of 
data from the fifth LIGO science run ^. The SNR 
threshold for SPA is 5.5, while for BCV-spin it is 8. With 
these thresholds, SPA and BCV-spin banks have approx- 
imately the same false alarm rates. 



B. EfHciencies 

We are now in a position to compare the efficiencies of 
SPA and BCV-spin banks. Given a large number of tar- 
get waveforms injected in the data, the efficiency is the 
ratio of the number of found injections to the total num- 
ber of injections made. For our purposes, an injection 
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Table II: Average number of triggers per time-shifts {{N)), and variances thereof {a = (iV^)^/^), for SPA with M = 0.95, 
and BCV-spin with stochastic placement and M = 0.8. Note that the 1-a intervals overlap, so that the false alarm rates are 
comparable. Next to these we list distances at which the efficiencies are 50%, 75%, and 90%. See [4^ for histograms of trigger 
numbers and efficiency plots. 



Bank 


(TV) 


a 


D^ox (Mpc) 


D75% (Mpc) 


Dgo% (Mpc) 


SPA 


97.3 


8.7 


40.1 


33.9 


15.9 


BCV-spin 


85.4 


8.4 


34.6 


17.5 


14.5 



is considered found if it had an SNR above the chosen 
threshold with at least one template in the bank, within 
a certain time interval around the time when the injec- 
tion was actually made. In the case of SPA, the width 
of this interval can be chosen to be 40 ms. BCV-spin 
templates, being phenomenological, turn out to have a 
larger timing inaccuracy, and an interval of 100 ms was 
found to be more appropriate. This had already been no- 
ticed in (4^ : presumably the larger timing uncertainties 
of BCV-spin are related to its unphysical phasing (es- 
sentially, missing PN terms) as it is predominantly the 
phasing which affects timing errors. 

It is important that efficiencies be compared for the 
same false alarm rate. And indeed, as we have just seen, 
SPA and BCV-spin have essentially the same background 
rates if the SNR thresholds are set at 5.5 and 8, respec- 
tively. 

We made 1124 injections distributed logarithmically 
in distance between 1 Mpc and 50 Mpc, with compo- 
nent masses randomly chosen between 1 Mq and 30 Mq, 
spin magnitudes |Si|/mf , z = 1, 2 between 0.7 and 1, and 
arbitrary directions for the initial spin vectors. The effi- 
ciency of SPA then came out to be 0.93, versus 0.89 for 
BCV-spin. These results have been summarized in Ap- 
pendix I of [4^ ; here we have provided a detailed account 
of how they were obtained. We refer to the latter paper 
for plots of efficiency against distance; see Table |ll] for 
the distances at which the efficiencies are 50%, 75%, and 
90%, both for SPA and stochastic BCV-spin. 

We find that despite the fact that BCV-spin banks 
have higher effective fitting factors with the target wave- 
forms than SPA banks, in a more realistic data-analysis 
comparison the two waveform families have similar abil- 
ities to detect simulated signals. The detection statis- 
tic for BCV-spin involves more degrees of freedom and 
the pipeline using BCV-spin waveforms is more sensitive 
to non-stationary noise transients in the data. Conse- 
quently, at the same false alarm rate the detection thresh- 
old of the BCV-spin bank is higher than the SPA bank, 
negating the effect of the improved effective fitting factor 
of the BCV-spin bank'^. Searches for spinning binaries us- 
ing the non-spinning bank therefore have approximately 
the same performance as even our improved BCV-spin 
bank. 



^ The problem had been anticipated in [4Cll |: here we have quanti- 
fied it using real data. 



VII. CONCLUSIONS 

Past searches for low mass compact binary inspiral 
events in LIGO data (with the exception of [431) have 
used waveforms which do not attempt to model the spin 
effects, despite the fact that astrophysical black holes 
may be spinning rapidly. In this paper we have con- 
structed template banks using the BCV-spin waveform 
proposed in [40[. Though phenomenological, these wave- 
forms seek to capture the spin-induced amplitude mod- 
ulation one expects to see in a physical signal, and have 
high effective fitting factors with PN waveforms that in- 
clude spin. We have improved on the search method 
of [4^ . in two ways: (i) we have constructed a bank us- 
ing the metric outlined in Ref . [4l| , which is much better 
suited to the template family, and (ii) we have explored 
two new placement algorithms (square-hexagonal and 
stochastic). We used spinning PN signals to study the 
effective fitting factors of three different banks: an SPA 
bank, a BCV-spin bank with square-hexagonal place- 
ment, and a BCV-spin bank with stochastic placement. 
We found that the two BCV-spin banks had a similar per- 
formance, but both did markedly better than SPA. How- 
ever, search performance should be judged by detection 
efficiency at a given false alarm rate. The search pipeline 
for low-mass compact binaries {2Mq < M < 35 M0) in 
data from the fifth LIGO science run used a non-spinning 
SPA bank with an SNR threshold of 5.5. We have demon- 
strated that to achieve a comparible false alarm rate with 
the currently available search pipelines using BCV-spin 
templates requires an SNR threshold of 8 and with this 
higher threshold, the detection efficiency of BCV-spin for 
spinning PN signals becomes similar to that of the non- 
spinning SPA pipeline. Our findings, presented at length 
here and summarized in Ref. [45| , were used to guide the 
decision not to repeat the analysis of Ref. [l^l with data 
from the fifth LIGO science run. 

In conclusion, the detection performance of the BCV- 
spin pipeline is similar to that of the non-spinning SPA 
pipeline. We note, however, that our comparison is not 
entirely fair, because the SPA pipeline implements the 
metric-based coincidence algorithm of Ref. [t^ which 
dramatically reduces the number of spurious coincident 
triggers. In principle such a technique could also be ap- 
plied for BCV-spin, but since the metric has essentially 
no dependence on (V'OjV'3) E^nd only a weak dependence 
on /3 it is unlikely that implementation of the metric- 
coincidence algorithm would improve the sensitivity of 
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the BCV-spin pipeline. This justifies the use of non- 
spinning SPA pipeUnes rather than BCV-spin pipeUnes 
in LIGO searches. Nevertheless, to search for spinning 
signals with non-spinning banks is still sub-optimal, and 
work is ongoing to improve the performance of searches 
for spinning signals using templates determined by physi- 
cal (rather than phenomcnological) parameters proposed 
in Ref. (btI . Ft^ . In the mean time wc recommend the 
continued use of non-spinning SPA banks in upcoming 
searches until more efficient template families designed 
to capture spin-modulated waveforms have been incor- 
porated into a pipeline. 



for many helpful discussions, and in particular Michele 
Vallisneri for a critical reading of the manuscript. CVDB, 
TC, IH, GJ, and BSS are supported by PPARC grant 
PP/FOO 1096/1. DB is supported by National Science 
Foimdation grant NSF-0847611. Hideyuki Tagoshi is 
supported by KAKENHI, Nos. 16540251 and 20540271. 
Hirotaka Takahashi is partially supported by the Uchida 
Energy Science Promotion Foundation, Sasagawa Scien- 
tific Research Grant and the JSPS Grant-in- Aid for Sci- 
entific Research No. 20540260. 



Acknowlegements 

We would like to thank the members of the LIGO Sci- 
entific Collaboration Compact Binary Coalescence group 



[1] B. Abbott et al. (LIGO Scientific Collaboration) (2007), 
0711.3041. 

[2] F. Acernese et al. (VIRGO), AIP Conf. Proc. 794, 307 
(2005). 

[3] B. Willke (LIGO Scientific Collaboration), Class. Quant. 

Grav. 24, S389 (2007). 
[4] B. Abbott et al. (2009), in preparation. 
[5] R. K. Kopparapu et al.. The Astrophysical Journal 675, 

1459 (2008), 0706.1283. 
[6] V. Kalogera et al., Astrophys. J. 601, L179 (2004), astro- 

ph/0312101. 

[7] R. O'Shaughnessy, C. Kim, V. Kalogera, and K. Belczyn- 
ski, Astrophys. J. 672, 479 (2008), astro-ph/0610076. 

[8] R. O'Shaughnessy, C. Kim, T. Frakgos, V. Kalogera, and 
K. Belczynski, Astrophys. J. 633, 1076 (2005), astro- 
ph/0504479. 

[9] L. A. Wainstein and V. D. Zubakov, Extraction of signals 
from noise (Prentice-Hall, Englewood CliflFs, NJ, 1962). 

[10] K. S. Thorne, in Three hundred years of gravitation, 
edited by S. W. Hawking and W. Israel (Cambridge Uni- 
versity Press, Cambridge, 1987), chap. 9, pp. 330-458. 

[11] L. Blanchet, T. Damour, and B. R. Iyer, Phys. Rev. D51, 
5360 (1995), gr-qc/9501029. 

[12] C. M. Will and A. G. Wiseman, Phys. Rev. D54, 4813 
(1996), gr-qc/9608012. 

[13] L. Blanchet, B. R. Iyer, C. M. Will, and A. G. Wiseman, 
Class. Quant. Grav. 13, 575 (1996), gr-qc/9602024. 

[14] L. Blanchet, Phys. Rev. D54, 1417 (1996), gr- 
qc/9603048. 

[15] L. Blanchet, T. Damour, B. R. Iyer, C. M. Will, and 
A. G. Wiseman, Phys. Rev. Lett. 74, 3515 (1995), gr- 
qc/9501027. 

[16] L. Blanchet, G. Faye, B. R. Iyer, and B. Joguet, Phys. 

Rev. D65, 061501 (2002), gr-qc/0105099. 
[17] L. Blanchet, T. Damour, G. Esposito-Farese, and B. R. 

Iyer, Phys. Rev. D71, 124004 (2005), gr-qc/0503044. 
[18] L. Blanchet, T. Damour, G. Esposito-Farese, and B. R. 

Iyer, Phys. Rev. Lett. 93, 091101 (2004), gr-qc/0406012. 
[19] K. G. Arun, L. Blanchet, B. R. Iyer, and M. S. S. 



Qusailah, Class. Quant. Grav. 21, 3771 (2004), gr- 
qc/0404085. 

[20] L. Blanchet, G. Faye, B. R. Iyer, and S. Sinha (2008), 
0802.1249. 

[21] L. E. Kidder, C. M. Will, and A. G. Wiseman, Phys. 

Rev. D 47, 4183 (1993), gr-qc/9211025. 
[22] L. E. Kidder, Phys. Rev. D 52, 821 (1995), gr- 

qc/9506022. 

[23] G. Faye, L. Blanchet, and A. Buonanno, Phys. Rev. D74, 

104033 (2006), gr-qc/0605139. 

[24] L. Blanchet, A. Buonanno, and G. Faye, Phys. Rev. D74, 

104034 (2006), gr-qc/0605140. 

[25] K. G. Arun, A. Buonanno, G. Faye, and E. Ochsner 

(2008), 0810.5336. 
[26] M. Hannam, S. Husa, U. Sperhake, B. Briigmann, 

and J. A. Gonzalez, Phys. Rev. D 77, 044020 (2008), 

0706.1305. 

[27] A. Buonanno, G. B. Cook, and F. Pretorius, Phys. Rev. 

D 75, 124018 (2007), gr-qc/0610122. 
[28] M. Boyle et al., Phys. Rev. D 78, 104020 (2008), 

0804.4184. 

[29] A. Buonanno et al., Phys. Rev. D 76, 104049 (2007), 
0706.3732. 

[30] A. Buonanno et al. (2009), 0902.0790. 
[31] B. J. Owen, Phys. Rev. D53, 6749 (1996), gr- 
qc/9511032. 

[32] B. J. Owen and B. S. Sathyaprakash, Phys. Rev. D 60, 

022002 (1999), gr-qc/9808076. 
[33] S. Babak, R. Balasubramanian, D. Churches, T. Coke- 

laer, and B. S. Sathyaprakash, Class. Quant. Grav. 23, 

5477 (2006), gr-qc/0604037. 
[34] B. Allen, W. G. Anderson, P. R. Brady, D. A. Brown, 

and J. D. E. Creighton (2005), gr-qc/0509116. 
[35] T. A. Apostolatos, C. Cutler, G. J. Sussman, and K. S. 

Thorne, Phys. Rev. D 49, 6274 (1994). 
[36] T. A. Apostolatos, Phys. Rev. D52, 605 (1995). 
[37] T. A. Apostolatos, Phys. Rev. D54, 2438 (1996). 
[38] V. Kalogera, Astrophys. J. 541, 319 (2000), astro- 

ph/9911417. 



11 



R. O'Shaughnessy, J. Kaplan, V. Kalogera, and [58 
K. Belczynski, Astrophys. J. 632, 1035 (2005), astro- 
ph/0503219. [59 
A. Buonanno, Y. Chen, and M. Vallisneri, Phys. Rev. 
D67, 104025 (2003), gr-qc/0211087. [60 

A. Buonanno, Y. Chen, Y. Pan, H. Tagoshi, and M. Val- 
hsneri, Phys. Rev. D72, 084027 (2005), gr-qc/0508064. [61 

B. Abbott et al. (LIGO Scientific Collaboration), Phys. 

Rev. D 78, 042002 (2008), 0712.2050. [62 
B. Abbott et al. (LIGO Scientific Collaboration), Phys. 
Rev. D73, 062001 (2006), gr-qc/0509129. [63 
B. Abbott et al. (LIGO Scientific Collaboration), Phys. 
Rev. D 77, 062002 (2008), 0704.3368. [64 
B. Abbott et al. (LIGO Scientific Collaboration) (2009), [65 
0901.0302. 

T. Cokelaer, Phys. Rev. D 76, 102004 (2007), 0706.4437. 

B. Allen, I. Harry, and B. S. Sathyaprakash (2009), in [66 

preparation. 

B. S. Sathyaprakash and S. V. Dhurandhar, Phys. Rev. [67 
D44, 3819 (1991). 

S. Droz, D. J. Knapp, E. Poisson, and B. J. Owen, Phys. [68 
Rev. D 59, 124016 (1999), gr-qc/9901076. 

B. Abbott et al. (LIGO Scientific Collaboration), Phys. [69 
Rev. D 69, 122001 (2004), gr-qc/0308069. 

B. Abbott et al. (LIGO Scientific Collaboration), Phys. [70 
Rev. D 72, 082001 (2005), gr-qc/0505041. 

B. Abbott et al. (LIGO Scientific Collaboration), Phys. [71 
Rev. D 72, 082002 (2005), gr-qc/0505042. 

B. Abbott et al. (LIGO Scientific Collaboration), Phys. [72 
Rev. D73, 102002 (2006), gr-qc/0512078. 
K. Belczynski, R. Taam, E. Rantsiou, and M. Van der 
Sluys (2007), astro-ph/0703131. [73 
G. B. Cook, S. L. Shapiro, and S. A. Teukolsky, Astro- 
phys. J. 424, 823 (1994). [74 
K. S. Thorne, Astrophys. J. 191, 507 (1974). 
L. Blanchet, Living Rev. Rel. 5, 3 (2002), gr-qc/0202016. 



T. Damour, B. R. Iyer, and B. S. Sathyaprakash, Phys. 
Rev. D63, 044023 (2001), gr-qc/0010009. 
L. Blanchet, B. R. Iyer, and B. Joguet, Phys. Rev. D65, 
064005 (2002), gr-qc/0105098. 

T. Damour, P. Jaranowski, and G. Schaefer, Phys. Rev. 

D 62, 084011 (2000), gr-qc/0005034. 

L. Blanchet, T. Damour, and G. Esposito-Farese, Phys. 

Rev. D69, 124007 (2004), gr-qc/0311052. 

L. Blanchet and B. R. Iyer, Phys. Rev. D71, 024004 

(2005), gr-qc/0409094. 

L. Blanchet, T. Damour, and B. R. Iyer, Class. Quant. 

Grav. 22, 155 (2005), gr-qc/0410021. 

H. Wahlquist, Gen. Rel. Grav. 19, 1101 (1987). 

K. G. Arun, B. R. Iyer, B. S. Sathyaprakash, and 

P. A. Sundararajan, Phys. Rev. D71, 084008 (2005), gr- 

qc/0411146. 

L. Lindblom, B. J. Owen, and D. A. Brown, Phys. Rev. 
D78, 124020 (2008), 0809.3844. 

Y. Pan, A. Buonanno, Y. Chen, and M. Vallisneri, Phys. 

Rev. D69, 104017 (2004), gr-qc/0310034. 

A. Buonanno, Y. Chen, and M. Vallisneri, Phys. Rev. 

D67, 024016 (2003), gr-qc/0205122. 

S. Babak, Class. Quant. Grav. 25, 195011 (2008), 

0801.4070. 

C. Messenger, R. Prix, and M. A. Papa (2008), 
0809.5223. 

LSC algorithm library, URL 

'http: //www. Isc- group. phys .uvnn. edu/lal| 
The LIGO Scientific Collaboration, Tech. Rep. 
LIGO-T070109-01, LIGO Project (2007), URL 
http : //www . ligo . caltech . edu/docs/T/T070109-01 . pdf] 
C. A. K. Robinson, B. S. Sathyaprakash, and A. S. Sen- 
gupta (2008), 0804.4816. 

A. Buonanno, Y. Chen, Y. Pan, and M. Vallisneri, Phys. 
Rev. D70, 104003 (2004), gr-qc/0405090. 



